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Preliminary results of a study of quenched chiral logarithms at /3 = 5.7 are presented. Four independent 
determinations of the quenched chiral log parameter 6 are obtained. Two of these are from estimates of the r;' 
mass, one from the residue of the hairpin diagram and the other from the topological susceptibility combined 
with the Witten-Veneziano formula. The other two determinations of S are from measurement of virtual i]' loop 
effects in vs. quark mass and in the chiral behavior of the pseudoscalar decay constant. All of our results are 
consistent with 5 — .080(15). The expected absence of quenched chiral logs in the axial-vector decay constant is 
also observed. 



1. Quenched Chiral Logs 

One of the most distinctive physical effects of 
hght quark loops in QCD is the screening of topo- 
logical charge, which is responsible for the large 
mass of the rj' meson. The absence of screening in 
the quenched approximation gives rise to singular 
chiral behavior (quenched chiral logs) [|l| arising 
from soft rj' loops. In full QCD rj' loops are fi- 
nite in the chiral limit, while in quenched QCD 
they produce logarithmic singularities which ex- 
ponentiate to give anomalous power behavior in 
the limit — s- 0. The first evidence of quenched 
chiral logs in the pion mass was reported last 
year|^,||. Here, and in a forthcoming paper [Q, 
we present the results of a more detailed study of 
this effect. 

In the context of the effective chiral Lagrangian 
for QCD, the ;7(3) x f7(3) chiral field U = 

exp iY^i=o't'i''^il f-^ exhibits singular chiral be- 
havior in the quenched theory arising from loga- 
rithmically divergent 77' — (j)Q loops: 



f/^exp[-((/)2)/2/2]c7 
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where U is finite in the chiral limit, and A is an 
upper cutoff on the 7/ loop integral. (We con- 



sider, for simplicity, the case of three equal mass 
light quarks.) The anomalous exponent S is de- 
termined by the ij' hairpin mass insertion uiq, 
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(Here is normalized to a phenomenological 
value of 95 MeV.) From the chiral behavior (|^) 
of the field U, we can infer the singularities ex- 
pected in the matrix elements of various quark 
bilinears. Matrix elements of operators which 
flip chirality, such as the pseudoscalar charge 
'ip^^ip (x U — U"^ should exhibit a chirally singu- 
lar factor (m^)~*. By contrast, operators which 
preserve chirality, such as the axial vector current 
tpj^-ff^ij; oc iU~'^d'^U + h.c are expected to be fi- 
nite in the chiral limit. From PCAC, it follows 
that 



777,^ oc TTlo^'' 



(3) 



In the data presented here, we verify all of these 
predictions by studying the pion mass and the 
vacuum-to-one-particle matrix elements 



(o|7/;7'^kb)) = fp 
(o|VV7''V'k(p)) = fAP'' 



(4) 
(5) 



*TaIk presented by H. Thacker 



The modified quenched approximation (MQA) 
[] provides a practical method for resolving the 



2 





16 




0.14 




0.12 




0.1 


C\J 
< 




CO 
CO 


0.08 


a 




E 




c 


0.06 


o 




"q. 






0.04 




0.02 



OL ^ ^ ^ ^ , ^ ^ 1 

5 10 15 20 25 30 35 40 
bare quark mass (latxl 0'*3) 

Figure 1. Chiral log effect in vs. niq. Solid 
line is a perturbative (quadratic) fit to the five 
largest masses (four are not shown). Dotted line 
is fit of all masses to Eq. (3). 



problem of exceptional configurations, and allows 
an accurate investigation of the chiral behavior 
of quenched QCD with Wilson-Dirac fermions. 
The pole-shifting prescription for constructing 
improved quark propagators is designed to re- 
move the displacement of real poles in the quark 
propagator, which is a lattice artifact, while re- 
taining the contribution of these poles to contin- 
uum physics. We have used the MQA procedure 
to study the masses and matrix elements of flavor 
singlet and octet pseudoscalar mesons in the chi- 
ral limit of quenched QCD. A complete discussion 
of these results will be presented in a forthcom- 
ing publication. Here we present results from 300 
gauge configurations sd, /3 = 5.7 on a 12'^ x 24 lat- 
tice, with clover-improved quarks {Csw = 1-57) at 
nine quark mass values covering a range of pion 
masses from .2387(53) to .5998(17). 

In addition to observing 77' loop effects, we also 
present two direct estimates of the chiral log pa- 
rameter 6, one from a calculation of the 77' hair- 
pin diagram, and the other from a calculation 
of the topological susceptibility combined with 
the Witten-Veneziano relation. The results for 
S agree well with each other and with the expo- 
nents extracted from fp/fA and m^, thus giv- 
ing four independent and consistent determina- 



tions of 6. All four results fall within one stan- 
dard deviation of S = .080(15). This is about 
a factor of two smaller than the result S = .17 
expected from the phenomenological values of 
Too « 850 MeV and U = 95 MeV. The agree- 
ment between the four determinations of S indi- 
cates that relations imposed by chiral symmetry 
and the Witten-Veneziano formula are approxi- 
mately valid, even at /? = 5.7. Possible reasons 
for the overall suppression of the exponent S com- 
pared to phenomenological expectations will be 
discussed in Ref . Q . 

2. The rf hairpin mass insertion 

Following Ref. ||] we calculate 7^ quark loops 
using quark propagators with a source given by 
unit color-spin vectors on all sites. The statisti- 
cal errors are dramatically improved by the MQA 
procedure, allowing a detailed study of the time- 
dependence even at the lightest masses. We de- 
termine the value of TOq with a one-parameter fit 
to the overall magnitude of the hairpin correla- 
tor, assuming a pure double-Goldstone pole form. 
In the chiral limit, we find toqa = .601(30), or 
Too = 709(35) MeV using a"! = 1.18 GeV. (For 
unimproved Csw = fermions, we get a much 
smaller value of 464(24) MeV.) Using Eq. (|), 
we obtain the values for the chiral log parame- 
ter S. Extrapolating to the chiral limit, this gives 
6 — .068(8) if we use the unrenormalized lattice 
value for /^r. If a tadpole improved renormaliza- 
tion factor is included, this becomes 6 = .095(8). 

3. Topological susceptibility 

The fermionic method for calculating the topo- 
logical susceptibilitylQ] can be implemented with 
the same allsource propagators used for the hair- 
pin calculation. For each configuration, we com- 
pute the integrated pseudoscalar charge Q5. The 
winding number v = —imqQ^ is then obtained 
and the ensemble average (i/^) is calculated. In 
the chiral limit, this gives Xt = (188 MeV)'^. 
(Here we have used = 1.18 GeV.) Using 
the Witten-Veneziano formula (with unrenormal- 
ized /tt) to obtain toq, the chiral log parameter 
S = .065(8) is found. 
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Figure 2. fp vs. in lattice units. Solid line is 
a linear fit to the four largest masses. 
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Figure 3. /a vs. in lattice units. Solid line is 
a linear fit to the four largest masses. 



4. Quenched chiral logs in the pion mass 

The pion masses are obtained for nine kappa 
values over a range of hopping parameters from 
.1400 to .1428, with Csw = 1-57, corresponding 
to quark masses from roughly the strange quark 
mass down to about four times the up and down 
quark average (i.e. a pion mass of « 270 MeV). 
A perturbative (linear+quadratic) fit of as a 
function of ruq works well for the five heaviest 
masses, up to k = .1420. The value of for the 
lighter quark masses deviates significantly below 
this perturbative fit (see Fig. 1) showing clear 
evidence of a quenched chiral log effect. Fitting 
to the formula gives an anomalous exponent 
S = .079(8). 

5. Pseudoscalar and axial-vector matrix el- 
ements 

By a combined fit of smeared-local pseu- 
doscalar and axial-vector propagators, we obtain 
values for the decay constants fp and /a defined 
in (^-(^). The behavior of these two constants 
as a function of pion mass squared is shown in 
Figs. 2 and 3. We see that there is a very signif- 
icant chiral log enhancement at light masses for 
the pseudoscalar constant, but the axial- vector 
shows no chiral log effect, just as theoretical ar- 
guments predicted. From the phenomenology of 



0(p'*) terms in the chiral Lagrangian, it can be 
argued that the perturbative slopes of fp and /a 
in full QCD should be approximately equal. In- 
deed, if a chiral log factor is removed from fp, 
the remaining slopes seen in our lattice data are 
equal within errors. To obtain an estimate of 6, 
we fit the ratio fp/ fA to a pure chiral log factor: 

^ = const. X (ml)-^ (6) 
fA 

This fit gives S = .080(7). 
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